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Context. The expansion of the Universe is accelerating, as testified by observations of supemovae of type la as a function of redshift. 
Explanations of this acceleration are of two types: modifications of Einstein gravity or new forms of energy, coined dark energy. 
Aims. The accelerated expansion is explained here by a combination of Dvali-Gabadadze-Porrati (DGP) model gravity and Chaplygin 
gas dark energy. Both models are characterized by a length scale L which may be the same. 

Methods. The continuity equation for the combined model is derived in flat geometry, and solved by numerical methods. 
Results. The solution is shown to have the expected properties: at very small scales (a <k L) the energy density behaves as pressureless 
dust, at very large scales (a » L) as a cosmological constant. The modifications to the DGP model and the Chaplygin gas model 
occur for values of a near L. 

Conclusions. The results show an increase in the present dark energy density relative to the plain DGP model. 
Key words, cosmology- dark energy 
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1. Introduction 

The demonstration by SNela that the Universe is undergoing an 
accelerated expansion has stimulated a vigorous search of mod- 
els to explain this unexpected fact. Since the dynamics of the 
Universe is conventionally described by the Friedmann equa- 
tions which follow from the Einstein equation in four dimen- 
sions, all departures from the KCDM concordance model must 
be due to some modifications of the Einstein equation. 

The left-hand-side of the Einstein equation encodes the ge- 
ometry of the Universe in the Einstein tensor G^, the right-hand- 
side encodes the energy content in the stress-energy tensor T MV . 
Thus modifications to G^y imply some alternative gravitation, 
modifications in 7^,, involve new forms of energy densities that 
have not been observed, and which therefore are called dark en- 
ergy. 

A well-studied model of modified gravity is the Dvali- 
Gabadadze-Porrati (DGP) braneworld model (Dvali & al. |2000l 
in which our 4-dimensional world is a FRW brane embedded in 
a 5-dimensional Minkowski bulk. On the 4-dimensional brane 
the action of gravity is proportional to M 2 t whereas in the bulk 
it is proportional to the corresponding quantity in 5 dimensions, 
Ml. The model is then characterized by a cross-over length scale 
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such that gravity is a 4-dimensional theory at scales a <K L where 
matter behaves as pressureless dust, but gravity "leaks out" into 
the bulk at scales a » L and matter approaches the behaviour of 
a cosmological constant. To explain the accelerated expansion 
which is of recent date (z ~ 0.5 or a as 2/3), L must be of the 
order of 1 . 

An equally interesting model introduces a fluid called 
Chaplygin gas (Kamenshchik & al. l200T1 Bilic & al. 120021 fol- 



lowing work in aerodynamics (Chaplygin |1904| l. This addition to 
T M y is intriguingly similar to the DGP model in the sense that it is 
also characterized by a cross-over length scale below which the 
gas behaves as pressureless dust, and above which it approaches 
the behaviour of a cosmological constant. This length scale is 
expected to be of the same order of magnitude as the L scale in 
the DGP model. 

Both the DGP model and Chaplygin gas have problems 
with fitting present cosmological data (briefly summarized by 
Copeland & al. 2006|l. In a comparison with supernova data the 
DGP model is slightly disfavored, at the level of 2cr (Rydbeck & 
al. 120071 1. In the Chaplygin gas model the Jeans instability of per- 
turbations behaves like CDM fluctuations in the dust-dominated 
stage (a <sc L), but disappears in the acceleration stage (a » L). 
The combined effect of suppression of perturbations and non- 
zero Jeans length leads to a strong ISW effect and thus of loss of 
power in CMB anisotropies (Amendola & al. 120031 Bento & al. 
120031 . 

This has led to generalizations which are less motivated than 
the original models. Below we shall discuss a model which com- 
bines both models, motivated by the similarities in their asymp- 
totic properties and in the cross-over scales. 



(1) 2. Continuity equations 



The Friedmann equation for Hubble expansion in the DGP 
model may be written (Deffayet & al. 1200 11 ) 



H 1 



H 2 
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(2) 
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where k = 8nG/3, and p is the total cosmic fluid energy den- 
sity. In the following we shall set k — corresponding to a flat 
geometry. A + sign in front of the root term causes accelerating 
expansion. 
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The Friedmann equation for the rate of change in H is sim- 
plest written 



2a 



-3kp 



(3) 



where p is the total cosmic fluid pressure. Pressureless dust has 
p = 0, in the KCDM model the pressure is p^ = -p\, in 
the Chaplygin gas model p 9 = -A/p^, where A is a constant. 
Differentiating Eq. |2]) and using it to eliminate the second time 
derivative in Eq. ^ one obtains the continuity equation (or en- 
ergy conservation equation). 

In the ordinary FRW universe, the continuity equation for 
Chaplygin gas is 



pyia) - 3H\p tp (a) - 



P<p(a) 



= 0. 



(4) 



In the DGP universe the continuity equation for a generic 
energy density p(a) with pressure p(a) is 

at 1 \ 3 

PW = -3 - I p(a) + p(a) + ^2 I + — lp(a) - p(a)]. (5) 

Since ordinary matter does not interact with Chaplygin gas, one 
can derive separate continuity equations for the energy densities 
p m and p v , respectively. Thus ordinary matter obeys Eq. d5]l (5) 
with p — 0, our combined model for Chaplygin gas in a DGP 
universe obeys 
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PM = -3- \p,(a) - — + —j + _ \p,(a) + — 



3. Solutions 

The continuity equation for Chaplygin gas integrates to 



(6) 



P<p(a) 



B 



(7) 



where B is an integration constant. Thus this models has two free 
parameters. Obviously its limiting behaviour is 



P<p(a) • 



„3 



(B a/6 /B\ I/6 
for a « I -J , p v {a) oc -p for a » 1 — 1 .(8) 



The identification of the cross-over limit here with L in the DGP 
model implies B/A m L 6 . However, this cannot be substituted 
into the combined model, because the latter may have a different 
integration constant. Moreover, Eq. (|6]l can only be integrated 
numerically, and then no value can be substituted into it. 

In the case of pressureless dust the integration of Eq. |5]l with 
p — gives 



Pm(a) 



-3 
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(9) 



where C is an integration constant. Otherwise there is only one 
free parameter, L. 

Obviously p m (a) oc a~ 3 for small values of a/L. The density 
must also be positive for all a, thus C must be non-negative. One 
can fix the value of C by some boundary condition at a » L 
without affecting the a-dependence of p,„(a) noticeably in the 
cosmologically interesting region a < 1. For a/L > 1, p m (a) 
approaches an aT x quintessence-like behaviour or de Sitter-like 
behaviour. Ultimately the exponential term in Eq. ([9]) begins to 
dominate, leading to a phantom-like solution. By choosing the 



integration constant C sufficiently small one can shift the de 
Sitter behavior and subsequent phantom behavior into the distant 
future without affecting the shape of p m (a) in the region a < 1. 

The continuity equation |6]l cannot be integrated analytically, 
but only numerically. Then a boundary condition for p^{a) needs 
to be specified. In addition the solution depends on the two pa- 
rameters A and L. The properties of the solution can best be stud- 
ied graphically. 

The function a 3 p ip (a) approaches a constant for a — > 0, thus 
the a-dependence of p ¥ ,(a) is gT 3 for all choices of A and L, in 
agreement with the functions (|7]) and d9]l. To permit a comparison 
of p^(fl) and p,„(a), the latter is normalized to unity at a — 0, and 
a 3 p (p (a) is scaled to unity with an arbitrary constant at a — 0.01 
(because the numerical function is singular at a = 0). 

For a > L the a-dependence of p ¥ ,(a) slows down, approach- 
ing a "cosmological" constant. At some later point the solution 
starts to grow exponentially in a phantom-like behavior, just as 
was the case for p„,(a) in Eq. ([9]). This cannot be cured exactly 
because of the numerical nature of the solution to Eq. ([6]), but 
the problem can be shifted into the distant future by a proper 
choice of boundary condition. The exact position and function 
value at the point of the boundary condition only affects the nor- 
malization of a 3 p ip (a), but not its shape in the cosmologically 
interesting region a < 1 . 

To appreciate the changes brought to the DGP model by 
Eq. (|6]l we shall study the ratio 



R(a, L) = 



p<p(a,L,A) 
p m (a, L, C) ' 



(10) 



Ideally, the boundary condition in Eq. Q should be chosen so 
that the magnitude of the exponential term hidden in the nu- 
merical solution p^(fl) would exactly compensate C in Eq. (|9j, 
in which case R(a,L) would approach a constant value when 
a » L. However, to meet such a condition is not possible with 
a numerical solution, thus one has to resort to represent the ef- 
fect of both integration constants by a normalization factor, here 
chosen to give /?(0.01, L) = 1. 

The interesting results of this model are the A- and L- de- 
pendences in the region of observed acceleration 0.5 < a < 1 (or 
1 > z > 0). We have found that A acts as a mere factor, that can 
be considered as included in the normalization condition for the 
ratio < fT0[ >. (For this reason we have not defined R to be a function 
of A.) 

Figure 1 gives an example of the L- dependence: we plot 
the ratio R(a,L) in the range 0.01 < a < 1 for L = 1.0 (thick 
solid line), L = 1.2 (thin solid line), and L — 1.4 (dashed line). 
The L-scale has been chosen as the double of the a-scale: L = 1 
corresponds to a = 0.5 or z = 1. It is clear from the Figure 
that p v (a) entails corrections to the DGP-model density p m (a) of 
the order of several percent in the region a m 0.5 - 1, where 
accelerated expansion has been observed. 

4. Conclusions 

We have studied a model of dark energy in which the geometry 
of the Universe is described by a brane in a 5-dimensional bulk: 
the DGP model, and where Chaplygin gas enters as a component 
in the stress-energy tensor. 

This combined model depends on two parameters, the DGP 
model cross-over length L and the Chaplygin gas parameter A, 
the latter entering only as a normalization parameter. In the re- 
gion a as 0.5 - 1 where accelerated expansion has been observed, 
the model density exceeds corrections to the DGP density by 
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Fig. 1. The ratio R(a, L) of the density p v {d) of the combined model 
to the DGP density p m (a) for L = 1.0 (thick solid line), L = 1.2 (thin 
solid line), and L = 1.4 (dashed line). L — 1 corresponds to a = 0.5. 



several percent, as can be seen in Fig. 1. Correspondingly, the 
density parameter £2^ increases at fixed Q m . 

In this letter no fit to SNela data has yet been done. CMB 
anisotropy data can only be used when power spectra have been 
derived for the 5-dimensional bulk space implied here. It re- 
mains to be studied whether the various shortcomings of the 
DGP model and the Chaplygin gas model then are overcome. 



References 

Amendola, L., Finelli, R, Burigana, C. & Carturan, D. 2003, JCAP 0307, 005 
Bento, M. D. C, Bertolami, O & Sem, A. A. 2003, Phys. Rev D, 67, 063003 
Bilic, N., Tupper, G. B. & Viollier, R. D. 2002, preprint arXiv: astro-ph/0207423 
Chaplygin, S. 1904, Sci. Mem. Moscow Univ. Math. Phys., 21,1 
Copeland, E. J., Sami, M. & Tsujikawa, S. 2006, IntJ.Mod.Phys. D15, 1753 
Deffayet, D. 2001, Phys. Lett. B, 502, 199; Deffayet, D, Dvali, G. R. & 

Gabadadze 2002, Phys. Rev. D, 65, 044023 
Dvali, G. R., Gabadadze & Porrati, M 2000, Phys. Lett. B, 485, 208 
Kamenshchik, A., Moschella, U. & Pasquier, V. 2001, Phys. Lett. B, 51 1, 265 
Rydbeck, S., Fairbairn, M. & Goobar, A. 2007, preprint arXiv: astro-ph/0701495 



